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In this paper a composite scattering model, suitable for 
explaining the behavior of measured scattering cross sections 
of the ocean surface, will be presented. Furthermore, uti- 
lizing this scattering model, the spectrums of the small 
gravity, gravity-capillary, waves will be predicted for 
NASA/MSC, 13.3 GHz Scatterometer data. 


INTRODUCTION 


From the viewpoint of radar scattering at high fre- 
quencies, the ocean scene is best described by the composite 
surface function £(x,y) given by 



5 = + + + + * CD 

where £ s is the swell, ^ is the sea, £ 2 is the 
gravity-capillary structure, and so on. The solution of 
the electromagnetic scattering from the sea requires the 
probability density function of the height and the correla- 
tion function of each component of £ , as well as the 
joint probability densities of the components of £ . The 
interaction of the surface wind with each component of £ 
must be formulated to study the wind dependence of the 
scattering cross section. However, a mathematical model 
of £(x,y) as a function of surface wind velocity is not 
available. 
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Most of the ocean studies in the past have been devoted 
to developing models for ocean surface wave forecasts. 

Toward this goal low frequency gravity wave spectrums for 
fully developed seas for various windspeeds and fetches have 
been studied (ref. 1). The dependence of significant wave 
height, ^ 1/2 > on surface windspeed has also been reported 
(ref. 2). The measurement of high frequency gravity-capillary 
waves has been reported in only a few instances, with the most 
recent investigation reported by Dobson (ref. 3). 

Due to the lack of complete mathematical description of 
£(x,y), the composite model for sea surface scattering will 
be studied using the NASA/MSC scattering cross section (a Q ) 
data. 

SCATTERING THEORIES AND COMPARISONS 

Numerous approaches have been advanced to explain scat- 
tering from rough surfaces. The three theories which have 
received attention and show promise of efficient interpreta- 
tion of experimental data are as follows: 

• The Kirchhoff method 

• The small perturbation theory for slightly rough 
surfaces 

• The composite scattering theory 

In the Kirchhoff method the field scattered by rough 
surface is formulated according to Huygen's principle and 
is given by the Stratton-Chu integral. The total field and 
its normal derivatives are determined by tangent plane 
approximation on th& surface. These requirements generally 
restrict this method to locally flat surfaces. 



The small perturbation method is valid for large values 
of the angles of incidence and exhibits meaningful polariza- 
tion dependence. It is useful in the low frequency limit' 
and can therefore be applied to a class of slightly rough 
surfaces when very low frequencies are used. A comparison 
of theoretical and experimental results over slightly rough 
water surface has been given by Wright (ref. 4). As pointed 
out in his paper, the measured average backscattering cross 
sections are in good agreement with the calculated values. 

The depolarized return from slightly rough ocean surface 
has been obtained by Valenzuela (ref. 5) by using Rice's 
small perturbation method. A slightly different approach 
using the small perturbation theory has also been used for 
slightly rough surfaces by Bass et al. (refs. 6 and 7). 

It is interesting to examine closely what parameters 
are needed in the formulation of the electromagnetic problem. 
If only the large-scale structure is considered, the tangent 
plane approximation can be used. With this theory o Q can 
be expressed as a function of tan $ Q , the root-mean-square 
slope of the surface, in the high frequency limit. Only a 
few measurements of the sea slope distributions have been 
reported. The most widely used of these are the optical 
measurements of Schooley (ref. 8) and Cox and Munk (ref. 9). 
Figure 1 shows the value of B Q as a function of windspeed. 
The curves C(l), C(2) have been calculated using the spec- 
trums of Kitaigorodskii and Pierson and Moskowitz (ref. 10). 
The lower value of C(l) and C(2), as compared to other curves 
in figure 1, is attributed to the fact that these curves are 
not direction dependent (upwind, downwind, etc.) but involve 
all facets of the sea surface. 
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For comparison, a value of tan $ Q = 0.27 was taken 
corresponding to a 20.5-knot upwind speed as given in 
figure 1. The value of o vv was calculated for the 
Gaussian height probability function. The calculated value 
of the scattering cross section (using the Kirchhoff method) 
given by 


(a ) 

v vv . 


sec 2 - e - iR.nCO) | 2 exp 
tan 2 $ Q 



( 2 ) 


for e = 55 + j30.25 is plotted in figure 2. In equation (2), 
0 is the incidence angle, and R i;l ( 0) is the Fresnel reflec- 
tion coefficient for normal incidence. Also shown in this 
figure are the NASA/MSC, 13.3-GHz, F4L8R1, 21-knot, forebeam 
data. The evident disagreement is attributed to the fact 
that equation (2) is a limiting solution, and only one com- 
ponent of the composite surface is considered. 


It has often been suggested that, near the normal direc- 
tion for backscattering cross sections, scattering of the 
optics type (Kirchhoff method) predominates. In other direc- 
tions, however, the slight roughness on top of the large- 
scale roughness constitutes the major source of scattering. 
The scattering cross sections are calculated from the 
scattered component of the field. In view of this, Wright 
(ref. 11) and Guinard and Daley (ref. 12) ignore the effect 
of large. structure to account for the scattering at higher 
backscattering angles. The procedure by Wright, Guinard 
and Daley parallels that of Rice (ref. 13), Barrick and 
Peake (ref. 14), and Valenzuela. For Rice’s method the 
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backscattering cross sections for a slightly rough surface 
using first order terms are given in reference 14. 




- 4tfk 4 cos 4 0 |a * | 2 W(p,q) 


In equation (3) 

k Q = wave number of the incident radar energy 
6 = the incident angle 

W(p ,q) - the roughness spectral density of the surface, and 
p,q are radian wave numbers 


e - 1 


cos 0 + Ve - sin 6 


a = ( £ ~ 1) [ (e - 1) sin 2 6 + el 
Je cos 0 + - sin 2 0 j 

e = the complex dielectric constant of the surface. 

The scattering cross sections can be computed from 
equation (3) for the exponential and Gaussian surface height 
correlation function. For the Gaussian correlated surface, 
the result is 


(a ) 

vv y 


4kjh 2 * 2 cos'* © | ot vv | ^ exp (-k 2 * 2 sin 2 0) 


2 n 2 .2 


where l is the surface height correlation distance, and 
2 . 9 
h is the surface-mean-square height. 
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Several ocean wave spectrums have been proposed over 
the last few years. However, no experimental spectrums of 
the fine structure (the capillary waves) have been reported 
for different wind velocities. In the equilibrium range the 
isotropic spectrum is of the form 


W(r) = 


12 , 2 

I p + q 


There is some uncertainty for the value of k , but 
Phillips (ref. 15) gives the following estimates: 

-3 

6x10 for the equilibrium range spectrum for gravity 
waves 

1.5x10 for the capillary wave spectrum 

Evaluating at wave numbers satisfying the Bragg scattering 
condition (ref. 11) equation (3) yields the limiting form 
of the cross sections as 




= 1 . 5tt x 10' 


a y6 | 2 cot 4 0 


In the case of the ocean, it is thought that the sea 
waves, the small gravity waves, and the gravity-capillary 
structure combined produce the scattering. The swells are 
assumed absent. To the first order of approximation, the 
composite scattering cross sections were calculated by adding 
the average incoherent scattering cross sections from the 
very rough surface (Kirchhoff method) to that of the slightly 
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rough surface. The mathematical proof of this is given in 
papers of Semenov (ref. 16) and Fuks (ref. 17). The com- 
parison of composite scattering cross section and NASA/MSC, 
13.3-GHz, F4L8R1, forebeam vv data is shown in figure 3. 

The theoretical curve is the summation of (a ) and fa 1 

vv l v vv's 

as given in equations (2) and (4) respectively. The value 
of the dielectric constant e is taken as 5S+j30.25 and, 
furthermore, tan B 0 = 0.27 . The result is encouraging. 
Comparisons such as shown in figure 3 made it obvious that 
the scattering by small gravity-capillary structure plays a 
significant role at higher angles. 

A comparison of equation (6) with 13.3GHz, NASA/MSC data 
in the range of angles 20° < 0 < 50° showed that the angular 
variation of the data was approximately the same as that given 
in the equation. It was therefore concluded that the direc- 
tional spectrums of the small gravity and gravity-capillary 
structure of the sea could be expressed as 

-k 3 

W(r) = kr J ( 7 ) 

The values of k and k 3 are wind dependent. After substi- 
tution of equation (7), and toward the goal of studying the 
change of spectrum as a function of wind velocity, an 
expression of the following form was used. 

%( 0 ) = k 1 W 1 2 |a vv | 2 (cos 0) 4 (cosec 0) 3 (8) 

In equation (8), is the wind velocity reduced to a 

19.5-meter anemometer height. After using equation (8) 
it was found that the value of (cos 0) 4 k 3 (in the range of 

14-7- 



angles 20° <01 50°) could be taken as 1 for most data. 
Consequently, the following simplified form of equation (8) 
was also used: 

k k 

a 0 (e) * k^ 2 1 Owl 2 (cot 6) 3 (9) 

By using algorithm 178 "direct search" from ACM communication 
(ref. 18) , a FORTRAN program was used to find the values of 
k lt k 2 , and k 3 . The program then searches for a minimum 
value. The value of dielectric constant was taken as 
e = 55+ j 30.25. 


Two typical sets of data will be analyzed here. One 
set consists of NASA/MSC Mission 119, 13.3 GHz data for 
vertical- transmit vertical-receive polarization combination. 
The upwind forebeam data, for F9L1R19 (flight 9, line 1 and 
run 19), F2L1R1 and F3L1R1, the corresponding surface wind 
velocities are 6 knots, 22.5 knots and 33 knots respectively, 
were processed. The values for the constants using equa- 
tion (8) for this set were as follows: 

k x = 0.026 

k 2 = v 1.324 
k 3 = 5.47 

In figure 4, the experimental and calculated data using 
preceding value of constants is presented. Similar results 
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were obtained for other sets of Mission 119 data. The same 
set of data using equation (9) gave the following values: 

k 1 = 0.043 

k 2 = 1.33 

k 3 = 5.00 

The second set of data was chosen from NASA/MSC Mis- 



sion 156 data. The data is the aft beam 13.3 GHz (vertical- 
transmit vertical-receive) data for upwind conditions. The 
F2L8R1, F5L4R5 and F6L4R2 data corresponding to 33 knots, 

15 knots and 3 knots of average wind speed respectively were 
processed. The values for constants using equation (8) were: 

k x = 0.0207 

\ 

k 2 = 1.1 

k 3 = 6.6 

Figure 5 shoAys the comparison of the experimental and cal- 
culated data using the preceding values of the constants. 

It should be pointed out here that the flight 6 data was 
for very calm conditions with extremely low surface wind 
velocity. 
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CONCLUSION 


It can be concluded that a theoretical composite model 
can explain the dependence of the scattering cross section 
on the angle of incidence 0 . The spectrum of the small 
scale structure is found wind dependent. In general, as 
the wind velocity increases, the magnitude of the spectrum 
for high spectrum frequencies increases. The value of the 
constant k 3 is also a function of the radar incident wave 
length. Two sets of data (F7L3R1, forebeam Mission 119 data) 
gathered under identical conditions for a surface wind of 
15.5 knots using equation (9) gave the following values for 


= 7.3 for 0.4 GHz , and 


4.7 for 13.3 GHz data . 


This dependence is expected since 0.4 GHz radar would be 
relatively insensitive to small gravity-capillary waves. 


14-10 




REFERENCES 

!• Ocean Wave Spectra. Proceedings of a conference 

sponsored by the U.S. Naval Oceanographic Office and 
the Division of Earth Sciences of the National Academy 
of Sciences, (National Research Council in 1961), 
Prentice Hall, N. J., 1963. 


m 





2. Moskowitz, L. : "Estimates of Power Spectrums for Fully 

Developed Seas for Wind Speeds of 20 to 40 Knots", 
Journal of Geophysical Research, Vol. 69, No. 24, 1964.. 

3. Dobson, E. B. : "Measurement of the Fine-Scale Structure 

of the Sea", Journal of Geophysical Research, Vol. 75, 
No. 15, May 1970. 

4. Wright, J. W. : "Backscattering From Capillary Waves 

With Application to Sea Clutter", IEEE Transactions on 
Antennas and Propagation, AP-14, No. 6, 1966. 

5. Valenzuela, G.: "Depolarization of EM Waves by Slightly 

Rough Surfaces", IEEE Transactions on Antennas and 
Propagation, AP-15, No. 4, 1967. 

6. Bass, F. G.; and Bocharov, V. G.: "On the Theory of 

Electromagnetic Wave Scattering From a Statistically 
Uneven Surface , Radio Engineering and Electronics, 

Vol. 3, No. 2, 1958. 

7. Bass, F. G., et al., "Very High Frequency Radiowave 
Scattering by a Disturbed Sea Surface", Parts I, II, 

IEEE Transactions on Antennas and Propagation, AP-16, 

No. 5, 1968. 

8. Schooley , A. H.: "A Simple Optical Method For Measuring 

the Statistical Distribution of Water Surface Slopes", 
Journal Opt. Soc. Amer., 44, 1954. 


14-11 








9. Cox, C.; and Munk, W. : Measurement of the Roughness of 

the Sea Surface From Photographs of Sun's Glitter", 
Journal Opt. Soo ., 44, 1954. 

10. Krishen, K. : "Saturation of o Q With Increasing Wind 

Velocity", Lockheed Electronics Company, HASD TR 

No. 649D. 21.007, Houston, Texas, 1969. 

11. Wright, J. W.: "A New Model for Sea Clutter", IEEE 

Transactions on Antenna and Propagation, AP-16, No. 2, 
1968. 

12. Guinard, N. W. ; and Daley, J. C.: "An Experimental 

Study of a Sea Clutter Model", Proc. IEEE Vol. 58, 

No. 4, 1970. 

13. Rice, S. 0.: "Reflection of Electromagnetic Waves From 

Slightly Rough Surfaces", Comm. Pure Appl. Math., Vol. 4, 
1951. 

14. Barrick, D. E.; and Peake, W. H.: "Scattering From 

Surfaces With Different Roughness Scales: Analysis 

and Interpretation", Research Report, Battelle Memorial 
Institute, Ohio, 1967. 

15. Phillips, 0. M.: The Dynamics of Upper Ocean , Cambridge 

University Press, London, 1966. 

16. Semenov, B.: "An Approximate Calculation of Scattering 

of Electromagnetic Waves From a Slightly Rough Surface", 
Radioteknika i Electronika (USSR), 11 (English Transla- 
tion) . 

17. Fuks, I.: "Contribution to the Theory of Radio Wave 

Scattering on the Perturbed Sea Surface", Izvestia 
Vyshikh Uchebnikh Zavedeniy Radiofizika, 5 (English 
Translation) . 


14-12 






18. Hooke, R. and T. Jeeves: "Direct Search Solution of 

Numerical and Statistical Problems", Journal Assoc. 
Computing Machinery t Vol. 8, No. 2, 1961. 


14-13 


Degrees 



Windspeed, knots 


Figure 1. - Values of 3 0 as a function of windspeed. 
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Figure 3. — Comparison of NASA/MSC, 13.3-GHz data with 

model theory. 
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